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Abstract 

This document provides a proof that the solutions to the convectively 
filtered Burgers equation, will converge to the entropy solution of the 
inviscid Burgers equation when certain restrictions are put on the initial 
conditions. It does so by first establishing convergence to a weak solution 
of the inviscid Burgers equation and then showing that the weak solution 
is the entropy solution. Then the results are extended to encompass more 
general initial conditions. 

1 Introduction 

Using a filtered velocity in fluid dynamics is not a new concept. Filtered 
velocities have been used in turbulence modeling in Large Eddy Simulation 
(LES) [11 [IIS], Lagrangian Averaged Navier-Stokes (LANS-a) [1 El El |7] , and 
Leray turbulence modeling 13 HO]- Specifically in the LANS-a and Leray 
approaches, a filtered velocity is used in the nonlinear term of the Navier-Stokes 
equations. A form of the compressible Euler equations with a filtered velocity 
has also been developed using the Lagrangian averaging [llj . In our earlier 
paper [T5], it was discussed that it should be possible to model both turbu- 
lence and shock formation using such a filtered velocity. This was motivated 
by realizing that turbulence and shocks are both consequences of the nonlinear 
term and its resulting cascade of energy into smaller scales. Thus it should be 
possible to capture both effects with proper small scale modeling. It has been 
seen that some turbulent behavior has been successfully modeled using a filtered 
velocity in the LANS-a and Leray approaches. This paper in conjunction with 
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our previous paper |12j , aims at showing that such a technique can successfully 
model shock formation. 

The investigation begins with the inviscid Burgers equation, 

ut + uux 0. (1) 

Burgers equation was chosen because it shares the same nonlinear term as the 
Euler and Navier-Stokes equations. Additionally it is a conservation law, like 
the Euler equations. It is known to form shocks, and has been well studied. 

It is well established that the inviscid Burgers equation forms discontinuities 
in finite time, determined by initial conditions [HI [14]. To deal with these 
discontinuities weak solutions are introduced. However, when weak solutions are 
introduced, solutions are no longer necessarily unique fl4l|15j . In order to choose 
the physically relevant solution, an entropy condition is applied, which one and 
only one weak solution satisfies. This physically relevant solution is referred to 
as the entropy solution. Lax, Oleinik, and Kruzkov have examined the entropy 
condition for conservation laws and expressed it using different techniques |14[ 
[T51 \W\ . Each of their entropy conditions can be used in different classes of 
conservation laws, but can all be applied to the inviscid Burgers equation with 
equivalent results 117|. This paper uses the Lax entropy condition, which is 
explained in section [2] 

Classically the inviscid Burgers equation is regularized by adding viscosity, 
resulting in the equation 

Ut + uu^^ vu^^. (2) 

This regularization has been proven to converge to the entropy solution of the 
inviscid Burgers equation as ^ [TU [T5J (TB] ■ 
This paper considers the equations 

Ut + uux = (3a) 
u — g°'*u (3b) 
u{x, 0) = uq{x) (3c) 

where 

5" = -.(-) (4) 
where g is a chosen filter. These equations replace the convective velocity of the 



inviscid Burgers equation with a filtered velocity. Thus, equations (3a I and (3b I 
are referred to as the convectively filtered Burgers equation (CFB). While it has 
been proven that the solutions to the CFB equations exist [12] , previously it has 
only been proven that the solutions for the Helmholtz filter converge to a weak 
solution of the inviscid Burgers equation with attempts to show numerically 
convergence to the entropy solution [15] . 

This paper proves that for a specific set of initial conditions that the solu- 
tions to the CFB equations will converge to the entropy solution of the inviscid 
Burgers equation. Specifically we will look at bell shaped, continuously differ- 
entiable initial conditions rigorously defined in Definition |4.1[ We then give 



2 



rationale and make a conjecture on how the CFB equations will converge to 
the entropy solution for any continuous initial conditions, and how to regain an 
entropy solution for discontinuous initial conditions. 

The following section reviews established facts about the inviscid Burgers 
equation and some of the recent work regarding the CFB equations. Section 
[3] proves that solutions to the CFB equations converge to a weak solution of 
the inviscid Burgers equation, and section [4] proves convergence to the entropy 
solution. Section [5] then extends the results of section |4] and conjectures that 
it can be extended further. Section [6] runs some numerical simulations and 
examines the results. All is the followed with concluding remarks. 

2 Background Information of Burgers Equation 
and the CFB equations 

Burgers equation has been thoroughly researched by many people over the years. 
This section provides a review of some of the previously established properties 
of the inviscid Burgers equation. Many of these will be used later on to establish 
new results about the CFB equations. This section will also list some of the 
previously established properties of the CFB equations, which are also crucial 
to the analysis found in the following sections. 



2.1 Method of Characteristics 

The inviscid Burgers equation lends itself well to examination with method 
of characteristics. From Whitham [13J, the inviscid Burgers equation can be 
broken into two ODE's 

MO^o (5) 

= (6) 

From this it is determined that along the characteristics 

^ = XQ + uoiOt (7) 

u{x) is constant. Thus characteristics travel at the speed equal to the value of 
u along those characteristics. This is, true, until characteristics cross, forming 
shocks. This is equivalent to seeing that the material derivative is zero [19]. 



2.2 Weak Solutions and Entropy Conditions 

Lax [14 addresses weak solutions and entropy solutions of conservation laws. 
From his work, a lot of information can be gained about the solutions to the 
inviscid Burgers equation. 
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The first thing we learn is that any weak solution to the inviscid Burgers 
equation must satisfy the integral form of the conservation law, or 



■h 



udx 




tl 



dt, 



(8) 



which must hold for any g and h and every time interval (ii, ^2)- A consequence 
of this are the Rankine-Hugoniot jump conditions. These dictate the speed at 
which any discontinuity can propagate. If s is the position of a shock then, 



Lax also establishes the existence and uniqueness of a weak solution to the 
inviscid Burgers equation which satisfies the so called entropy condition 



where s is the location of a discontinuity. Thus the only discontinuities that axe 
allowed to exist in this "entropy solution" are decreasing jumps. 

Lax also states that for solutions satisfying the entropy condition, "every 
point can be connected by a backward drawn characteristic to a point on the 
initial line." Thus any value of the entropy solution, u{x,t), can be traced back 
to the initial conditions. For discontinuous initial conditions, points traced back 
to the point of discontinuity can take on values between the left and right limits 
of the discontinuity as is shown in section |5.1| For continuous initial conditions 
the entropy solution can be written as u(x,t) = uo((j)(x,t)), where <j){x^t) is an 
increasing function of x for any time, and 4>{x, 0) = x. 

Here we will define what will be referred to in this paper as a reparameteri- 
zation of a function. 

Definition 2.1 // <j){x,t) is an increasing Junction of x for any time, and 
(j){x,Q) — X, the Junction J{<f){x,t)) will he called a reparameterization oj the 
Junction J. 

It is clear that at any time t, a reparameterization of the function / cannot 
obtain values that are not obtained by /. It can, however, lack values that are 
found in /, as it was not dictated that (j) be onto for all time. Looking back 
to the previous paragraph we can see that for continuous initial conditions, the 
entropy solution to the inviscid Burgers equation will be a reparameterization 
of the initial conditions. 

2.3 Properties of the CFB equations 

From previous work by our group [20j |2T1 [22j [12] the following theorem is es- 
tablished. It is presented here in its one dimensional form. 



(9) 



(10) 
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Theorem 2.2 Let g{x) S Vri'i(M) and uo{x) G C^(K), then there exists a 
unique global solution u{x,t) G C"'^(M,M) to the initial value problem (11). 



ut + uux — (US') 
u — g * u (lib) 
m(x, 0) = uq{x) (11c) 

A sketch of the theorem is as follows. Examine the equations using method 
of characteristics. Due to the nature of the equations the infinity norm of u will 
be bounded for all time. By Young's inequality ||u^||oo can thus be bounded 
for all time. The characteristics of the equations will not cross if their Jacobian 
remains nonzero. The rate of change of the Jacobian can be directly related to 
Ux by 

^^J = UxJ. (12) 

Since ||u2;||oo remains bounded, the Jacobian will remain nonzero, the charac- 
teristics will not cross, and a unique solution will exist for any finite time. 

In the course of proving the theorem, it was established that the solution 
take the form u(x,t) = ^□('/'(s;, i)) where 4>{x,t) is a continuous, invertible, and 
increasing function of x for any time, and (j){x,0) = x. Thus the solution is a 
reparameterization of its initial conditions. 



3 Weak Solution 

Regularizations of conservation laws do not necessarily have to converge to weak 
solutions to those conservation law. Take for example the KdV equations, 

Ut + uux ^ -euxxx- (13) 

This regularizes the inviscid Burgers equation in the sense that solutions are 
now continuous, however, many oscillations form as e ^ 0, requiring a weak 
limit for convergence [531 This limit is not a weak solution of the inviscid 
Burgers equation and thus definitely not the entropy solution. 

Thus the first step to proving convergence to the entropy solution is to prove 
convergence to a weak solution. The following subsections prove this by showing 
that a subsequence of the solutions to the CFB equations must converge to a 
function in L^^^. It is then shown that this function is, in fact, a weak solution 
to the inviscid Burgers equation. 

3.1 Convergence of Solutions 

In this section we show that the solutions of the CFB equations (u") converge 
to a function u. This subsection mirrors work done by Bhat and Fetecau |18j . 
We begin by claiming the following properties of the solutions u". 
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Lemma 3.1 The solutions to the initial value problem ^ have the following 
properties. 

\\u"{;t)\\L^ = 1 1?/" (-,0)1^00 = WuoWl^ = (14) 

TViu''{;t)) = TV{u^{-, 0)) = TV{uo{-)) = ^2, (15) 
\u°'{x,t) - u'^{x,s)\dx < Aslt - s|, (16) 

where A\, A2, and A^ are independent of a and TV{f{-)) can be defined for a 
smooth function f as 

TV{f{-))= [ \f'{x)\dx. (17) 



Proof Property 14 is verified by the existence proof in earlier papers [211 1221 [12| . 



that ||u"(-,i)||i=. = ||u"(-,0)|Uoe. 

To verify property |15[ take the derivative of 3a multiply by sign{u^) and 
integrate over the real line to obtain 



\ux\dx+ sgn{ux){uux)xdx = 0. (18) 



d_ 
dt 

Break the second term into intervals where sign(u^) remains constant, and 
u are continuous due to previous existence theorems, so at the locations that 
sign{ux) switches signs, the value of will be 0. Thus the second term is zero 
and we obtain the result 

\\u.{,-Ml^ = ll".(-,0)|Ui, (19) 

and thus Property [T5| is established. 

Property |16| can be proved by the following estimate: 



\u°'{x,t) -u°'{x,s)\dx < 



\u?\dtdx 



< AiA2\t-s 



rldt 



From Bressan [55] and Serre ^\ we know that properties (14 1, (15 1, and 
( 16 ) are enough to guarantee that a subsequence of u" converges to a function u 
in iiQp. Furthermore, u shares the same infinity norm bound as that established 
in (14 1, the same total variation bound as that in (151. 
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3.2 Convergence to a WeEik Solution 



To begin wc look at a specific subset of filters. The filters we examine are the 
functions whose Fourier transforms can be written as 



m = 



1 



Noting that gii = u we can see that 



with n < oo,Cj >0 C„^ 0. 



and 



02^ 



We will refer to a filter of this form as satisfying condition A. This class of filters 
includes the Helmholtz filter, which has been of previous interest in turbulence 
modeling. 

Clearly g{x) and its derivatives up to (yi(2"-2) (^2;) are well defined and bounded 
^ i+Y^" — c~fc2j absolutely integrable. 

If u and its derivative are absolutely integrable, then for a g satisfying 
condition A, the convolution 



tU = 



7 — 1 ^ j; 

is well defined. Furthermore, by Young's inequality 

1 











< 

00 


dxi-^^ 



111 — 



a-? 



,0-1) 



Thus there exists a constant A4 such that 



dx^' 



< -^A4 for j < 2n - 1 



This criteria is used in the following lemma. 

Lemma 3.2 Let it" be a sequence of functions that satisfy the following condi- 
tions. 



j K\dx,j \u'^\dx < A2 



dx^i 



< —A4 



for j <2n — 1 



(20a) 
(20b) 

(20c) 
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Let f S C°° be compactly supported on M. Then as a — > the quantity 



a'"/ (21) 



limits to 0. 



Proof For convenience the u°' shall be denoted u. Integrate Equation (21) by 
parts to obtain 



a^" y u^fdx = a^" j £J^(^^./) dx. 
Use product rule to expand (J^Wa;/) 



(22) 



(23) 



Take the absolute value, separate the last two terms of the binomial expan- 
sion, and apply the triangle inequality. 



< |a2"/u(")il("+i)/da;| 

2„ J ^{n) Y^n^2 ^(1+^) j(n-») 



Begin by bounding the third term, 
2« j ^in) ^V^')^z(i+0/("-)da; 



n-2 



i=0 



<a2"^ ;)ii^^"^iuii«^^+'^iuii/^ '111 



(24) 
(25) 
(26) 



f(n-i)| 



3rd term 



By applying the bound on ||u 



n^2 , s 

3rd term < ^ (^"j a"-'-iA4| l/^""^) | |i, 



(27) 



(28) 



which limits to as a — s- 0. 
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Next deal with the second term. 



2nd term 



nn — 1 



a n u 



n-l 

E 



9x' 



1 



n — 1 



(29) 



(30) 



(31) 



Again, apply the bound on ||u'-*-'||oo to get 



2nd term < n a" 

i=0 



n — 1 
i 



^4||^iW||l||/("-^^||c 



(32) 



(33) 



Since / and all its derivatives are bounded and < A2 the second term 

also limits to zero. 

Now for the first term . 



(34) 
(35) 



This differs from the second term only by a constant, so it must limit to 
as a ^ 0. 

Thus we obtain the result 



lim a^" 

a— >0 



00 / Q2n 



dx^ 



u ] Uxf dx — 0. 



(36) 



The last piece needed is taken from Duoandikoetxea [25]. The following 
lemma is a restatement of Duoandikoetxea' Theorem 2.1 from page 25. 

Lemma 3.3 Let g be an integrable function on M such that J g = I. Define 
9"=^9{^)- Then 

lim||5"*/-/||p = 
if f Cz U\l < p < 00 and uniformly ( i. e. when p = 00) if f £ Co(K) . 



With lemmas |3.3| and |3.2| we can now prove the following theorem regarding 
convergence to weak solutions. 
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Theorem 3.4 For any g satisfying condition A, the solutions u" to the CFB 
equations converge to a weak solution of the inviscid Burgers equation. 

Proof It was already shown that u" converges to a function u. To show this is 
a weak solution of the inviscid Burgers equation, we need to prove that for any 
test function / e C°° that has compact support on M x [0,r] that 

uft + -u'f^dxdt^O. (37) 



Begin by rewriting Equation 3a 

1 



as 



^(S")^) =(u^-<)^2". (38) 



Multiply by the test function / and integrate over M x [0,r] 



<f+{liu")^] fdxdt^ I I {W^ - f dxdt (39) 



Integrate by parts. 



' ""/t + d^dt = ( I {u'^- u'')u"f, dxdt 



2' 



{u" -u")u'^fdxdt (40) 



Taking the limit as a ^ of the left side, you get 

uft + fx dxdt 

Clearly then if the right hand side limits to zero, we have u is a weak solution 
to Burgers equation. 

Begin with the first term on the right hand side of Equation |40] The first 
term can be shown to limit to zero by noting that ||u"||oo has a uniform bound 
of Ai, and that since / G G°° with compact support, there exists an F € M"^ 
such that 1 1 /I loo < F and ||/a;||oo < F. Additionally let / be supported on the 
compact set fi. This leads to the bound 

^ y"(u"-u")u"/,dx(ii <FAiT||M"-u"||ii(n). (41) 

Take the limit of ||u" — ^"1^1(0). Break apart the norm with the triangle 
inequality to get 

lim - *.g"|| < Yan\\u'^ ~u\\ + \\u-u*g"\\ + \\g" *{u-u°')\\ 

a— »0 a^O 
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where the norms are all || • ||Li(a). The first and third term limit to zero as u° 



converges to u in L\^^. The second term limits to zero by lemma 



3.3 



Now deal with the second term from equation (40 1. Since g satisfies condition 



the second term can be rewritten as 

i-T 







By Lemma |3.2| every term in the sum limits to zero. Hence the sum limits 
to zero. 



Therefore the limit as a ^ of Equation ( 40 1 becomes 

/ '^ft + -u^fxdxdt = 0. 

"'K ^ 



(42) 



proving u is a weak solution of the inviscid Burgers equation. 



4 Convergence to the Entropy Solution 

In this section we will first examine some of the properties of non-entropic solu- 
tions, that is solutions that are a weak solution to the inviscid Burgers Equation, 
but do not satisfy the entropy condition. By examining these properties, it will 
be shown that the solutions to the CFB equation lack certain properties found 
in all non-entropic solution. Thus it will be shown that the solutions to the CFB 
equations converge to the entropy solution of the inviscid Burgers Equation. 

This examination will be limited to a class of initial conditions. Specifically, 
we intend to examine initial conditions that are continuously differentiable, and 
are bell shaped, i.e. have an interval where the functions are increasing, followed 
by an interval where the functions are decreasing. Functions that satisfy this 
condition will be referred to as satisfying condition B. It is for these functions 
as initial conditions that we will prove convergence to the entropy solution. 

Definition 4.1 Let u{x) £ C^(M) and Ux > over {—oo,p) and Ux < over 
{p, —oo) for some p. Additionally let u(x) have finite limits as x ±oo Then 
u{x) is said to have satisfied condition B. 

4.1 Non-entropic Weak Solutions 

There are three classic types of entropy violating weak solutions to the inviscid 
Burgers equation. This subsection shows examples of each type. The first 
is when you start with an increasing shock in the initial conditions and then 
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that shock remains, propagating at the speed dictated by the Rankine-Hugoniot 
jump conditions. An example of this is 



if a; < ^t, 

1 if U< x, 



u{x, t) 

taken from Lax |14j and is illustrated in figure [l] 



(a) 



(b) 



Figure 1: Here a discontinuity is introduced in the initial conditions and remains. 
The shock must travel at the speed dictated by the Rankine-Hugoniot jump 
conditions to be a weak solution. 



The second case is when a shock already exists and then splits into multiple 
shocks, one of which is an entropy violating shock. All the shocks move with 
the speed dictated by the Rankine-Hugoniot conditions. For a > 1 the following 
is a weak solution to Burgers Equation. This example was taken from Oleinik 
[T5] and is illustrated in figure [2] 

1 if a; < ^t, 
-a if < x < 
a if < a; < 
-1 if < X 

Another example is spontaneous shock formation with shocks forming out 
of a continuous interval. For a > the following is a weak solution to Burgers 
Equation. This example was taken from Serre [27] and is illustrated in figure [s] 

if a; < ^t, 
-a if < a; < 
a if < a; < |i 
if f i < X 

In the next section it is shown that these three cases exemplify the only type 
of entropy violating behavior possible. 





12 



0.B - 

U 0- 

-0.5 - 
-1 - 



(a) (b) 

Figure 2: A shock can split into multiple shocks and still remain a weak solu- 
tion, (a) The initial conditions, (b) The solution after the shock splitting has 
occurred. 



(a) (b) 

Figure 3: A shock can form from a continuous interval, (a) The initial condi- 
tions, (b) The solution after spontaneous shock formation has occurred. 



4.2 Decreasing Slope along characteristics 

By examining the inviscid Burgers Equation, it is possible to see that a non- 
entropic solution cannot form through the steepening of the solution. With this 
information we can then limit the ways a non-entropic solution can come into 
being. Begin with the inviscid Burgers Equation 

Ut + uuj. ~ 0. (43) 

In section [2] it was seen that along the characteristics 

^^xo + uoiOt (44) 
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that the value of u remains constant. This is true until characteristics crossed, 
at which point a shock is formed. 

Here a similar approach is taken, but on the derivative of the inviscid Burgers 
Equation. Differentiate the inviscid Burgers Equation to get 

^{ux) + u-^{ux) ^ -(u^)'^ (45) 
at ax 

Now if you examine this equation you find that along the same characteristics 
^ — xq + uo{£,)t that the quantity Ux is governed by 

= -{uxf. (46) 

at 

Thus for piecewise differentiable solutions, is always decreasing along charac- 
teristics and an increasing shock cannot form from the steepening of the solution. 

Now consider a solution that begins with initial conditions satisfying con- 
dition B. That solution is a continuously differentiable solution to the inviscid 
Burgers equation and is thus an entropy solution. It will remain an entropy 
solution until an increasing jump is formed. An entropy solution for initial con- 
ditions satisfying condition B will be piecewise continuous and thus from above 
will not steepen into an increasing shock. From this we conclude that an in- 
creasing shock can only occur if it exists in the initial conditions or must form 
instantaneously as it cannot form from the steepening of the solution. It can 
either form at existing points of discontinuity or form at points of continuity, 
which this paper refers to as shock splitting and spontaneous shock formation 
respectively. 



4.3 Entropy violating solutions are not reparameteriza- 
tions of initial conditions 

In section |2] it was established that the entropy solution of the inviscid Burgers 
equation is a reparameterization of initial conditions when the initial conditions 
are continuous. This subsection shows that a non-entropic solution cannot be 
both a weak solution and a reparameterization of initial conditions satisfying 
condition B. 

We first begin be examining some consequences of being both a weak solution 
and a reparameterization of initial conditions satisfying condition B. Then we 
assume that there is a non-entropic solution that is both a weak solution and a 
reparameterization and show that this is a contradiction. 

If a function is a reparameterization of initial conditions satisfying condition 
B, it is easy to see that the reparameterization will have one interval where it 
is increasing followed by an interval where it is decreasing. It is also clearly 
bounded. However, it need not be continuous. As a direct consequence of the 
Monotone Convergence Theorem for sequences, every point on the reparame- 
terization will have a well defined left and right sided limit. Since the left and 
right sided limits are well defined, the only type of discontinuity allowed is a 
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jump discontinuity. If a more rigorous explanation is desired we refer the reader 
to Section 5.7 in Davidson and Donsig [35]. 

Additionally any function satisfying condition B will have bounded variation. 
Thus any function that is a reparameterization will have variation bounded by 
the original function's variation. Thus if a solution is a reparameterization of 
initial conditions initial conditions satisfying condition B then it is of bounded 
variation. 

From Theorem 1.8.1 on page 21 and page 52 in Dafermos [30] we know that a 
function u that is of class BVioc and is a weak solution will satisfy the Rankine- 
Hugoniot jump conditions at every jump discontinuity. This means that if x is 
the location of a discontinuity then 

dt^ 2 ' ^ ' 

Thus if the solution is a weak solution and a reparameterization of initial con- 
ditions satisfying condition B, all of its discontinuities must be jump disconti- 
nuities satisfying the Rankine-Hugoniot jump conditions. 

To show that a function is not a reparameterization of initial conditions 
satisfying condition B, it is sufficient to find three points xi < X2 < such that 
u{xi) > u{x2) and u{x2) < u^x^). Essentially a function satisfying condition 
B is bell shaped and finding these three points finds an upsidedown bell, which 
cannot happen in a reparameterization. This is precisely the method used to 
show that a non-entropic solution cannot be both a reparameterization of the 
initial conditions. 

Since we are considering only initial conditions satisfying condition B, we 
are only beginning with continuous initial conditions. Thus from section [4?2| the 
only possibility of having a non-entropic solution is through spontaneous shock 
formation or by shock splitting. It will be shown that if either of these occur, 
then the non-entropic solution fails to be a reparameterization of the initial 
conditions. 

The following lemma is used later on when dealing with spontaneous shock 
formation and shock splitting. Because a non-entropic solution must still be a 
weak solution, spontaneous shock formation and shock splitting must behave is 
certain ways. The inviscid Burgers equation can be considered as a conservation 
law of wavemass, J u. Lemma 



4.2 



uses this fact to place restrictions on how 
spontaneous shock formation and shock splitting can occur. 

Lemma |4.2| addresses the area between the leftmost and rightmost shock, 
when a spontaneous shock formation or shock splitting occur. Essentially it 
says that if the area between the shocks has a higher value than the value on 
the outside of the shocks, then wavemass has been created, and it is no longer 
a weak solution to the inviscid Burgers equation. Figure |4] shows illustrations 
of this. 

The lemma proves that if the area between the leftmost and rightmost shock 
has values greater than those on its borders, then u(a;, t) cannot be a weak 
solution of the inviscid Burgers equation and a reparameterization of initial 
conditions. This is a proof by contradiction, so we assume that u{x,t) is a 
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1 

Extia Wavemas5 Is 

J 



(a) (b) 

Figure 4: Wavemass is created in shock splitting or spontaneous shock formation 
if the middle area is greater than the outer area's values at the shocks, (a) The 
conditions of u{x,t) before shock splitting. (b)If a shock splits and the middle 
area has higher values than its surroundings, extra wavemass has been created. 



Here xi Xr are the left most and right most shocks as used in lemma4.2 
The extra wavemass is indicated. 



weak solution and a reparameterization of initial conditions satisfying condition 
B, which places several constraints on it. Such a weak solution to the inviscid 
Burgers equation must satisfy the Rankine-Hugoniot jump conditions, Equation 

Additionally, weak solutions must satisfy the integral form of the conserva- 
tion law 



idx 



or if g and h are moving boundaries 

h(t) *2 ^ 2 ''(*) 



u dx 



9(t) 



tl 



9(t) 



dt 



dt. 



h{t) u{h{t),t) 



(48) 



d 

~dV 



g{t)]u{g{t),t) dt. 



The second definition. Equation |49[ is used in the following lemma. 
Lemma 4.2 Assume that u{x, t) takes the form 

{a{x,t) ifx<xi{t), 
b{x,t) if Xi{t) < X < Xr{t) 
c{x,t) ifXr{t)<x, 



(49) 



(50) 



where X/(^) ^^'^ Xr(t) are locations of discontinuities and X;(ti) = Xr(^i) = x* ■ 
At time ti let a(x*^,ti) > c{x*~^,ti). If for some period of time after ti and all 
X G iXhXr), b{x,t) > a{xi{t)^ ,t) and h{x,t) > c{xr{'t)^ ,i) , then u{x,t) cannot 
he a weak solution of the inviscid Burgers equation. 
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Proof Begin by assuming that u{x, t) is a weak solution of the inviscid Burgers 
equation and thus must satisfy Equation [49| for any g{t) and h{t). We will start 
by considering the left hand side of Equation|49]with selected moving boundaries 
and show that it is strictly greater than the right hand side, proving that u{x, t) 
cannot be a weak solution by contradiction. The moving boundaries will be 
defined by the positions of the leftmost and rightmost shock. 

With the moving boundaries established, begin with the left hand side of 
Equation (49 1 for the given boundaries. By putting a bound on the integrand, 
we transform the spatial integral into a temporal integral. 

u dx 



Xr(t) * = 



Xl{t) 



> min b{x,t2){xr-Xi) (51) 
ti 



LHS 

d d 

min HxM) / -KiXr - TT^Xi dt. (52) 



Now manipulate the equation to begin resembling the right hand side of 
Equation |49) 

f d d \ 

LHS > / min b{x,t2) -^Xr - tt^Xi dt 
Jti ^6(xiX.) \dt dt J 

*V . w ^ . .\ f d d 

laiii b{x,t2)-c\Xr,t)\\^Xr-^,Xi 

u Veixixr) J \c>t at 

d 

+ i^iXht) - c{Xr,t)) -7^X1 

+ c{Xr,t)—Xr - a{Xht)g^Xi dt 
Substitute in the speed of xi dictated by the Rankine-Hugoniot jump conditions. 

( b{xht) + a{Xhty 



+ (a(x;,i) - c{xr,t)) 



I \ 9 , . d , 

+ c{Xr,t)—Xr - a[xut)—xi dt 

^ . ^ . ^^ / 5 d 

mm 6(x, t2) ~ c[Xr, t) — Xr - -^.Xi 
xe(xiXT) J \ot dt 



+ (a(x;,i) - c{xr,t)) 



L 

f b{xi,t) - c{Xr,t) 
\ 2 



M 



+ {a{xi,t)-c{xr,t)) (^ "^'^'•' + ) + cixr, t)^^xr - a{xi,t)^^Xi 
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Consider term L. The value of b{x,t) for all x and some period of time after 
ti was designated to be higher than c{xr,t)- Additionally, for at least a short 
period of time §^Xr > §iXh otherwise the interval {xi , Xr) cannot have a nonzero 
measure. Thus for values t2 close to ti, term L is strictly positive. 

Now consider term M. Again the value of b{x, t) for all x and some period of 
time after ti was designated to be higher than c{xr,t). It was designated that 
at time ti, a{x*~,ti) > c{x*~^,ti). If a{x*~,ti) > c{x*'^,ti), then for values 
t2 close to ti term M is strictly positive. If a{x*~,ti) = c{x*^,ti), then by 
choosing t2 close to ti, term M can be made arbitrarily small. 

As t2 approaches ti, term L is approaching a strictly positive number, and 
term M is approaching a non-negative number. Thus it is possible to choose a 
t2 where J^^ L + Mdt > 0. Using this we see that 

rrr^^ a{Xutf-c{Xr,tf , d . ,^ ^ 

LHS>J l-cixr^tj—Xr-aixhtj—Xidt. (53) 

The right hand side of Equation [53] is the right hand side of Equation [49] with 
our chosen boundaries. Since with our moving boundaries the left hand side 



of Equation 49 is strictly greater than the right hand side, u{x, t) cannot be a 
weak solution. 

This result is now used to show that if there is spontaneous shock formation 
or shock splitting that u{x, t) cannot be both a weak solution and a reparame- 
terization of initial conditions. 

4.3.1 Spontaneous shock formation 

Assume that u(x, t) is the entropy solution to the inviscid Burgers equation 
up to time ti where an increasing shock is formed spontaneously at point x* . 
For such a discontinuity to form at least one other discontinuity must form in 
response. Thus we say after time ti, 

if x < xiit), 

u{x, t) = { b{x, t) if xi{t) <x< Xr{t) (54) 

if Xr{t) < X, 

where a(x, t), b{x, t), and c{x, t) are weak solutions to the inviscid Burgers equa- 
tion and xi{t) and Xr{i) give the position of the leftmost and rightmost discon- 
tinuities formed during the shock splitting. There may be more than two shocks 




formed as seen in section (4.1), but we just need to examine the leftmost and 
right most. 

For u{x, t) to be a weak solution we note that several things must be true. 
The speed of xK^) ™d Xr{t) are dictated by the Rankine-Hugoniot jump con- 
ditions to be 

d _ a{xT)+Kxt) d _ bjxr) + c(x+) , . 

j^Xi- ^ j^Xr- ^ . (55) 
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For there to be spontaneous shock forming, there must be some interval (ti, t2), 
where -^Xr > ^Xi- Thus for some interval {ti,t2), if aixT) — then 
b{xiit)~^,t) < b{xr{t)^,t)- Assume that t e (^1,^2) for the remainder of the 
subsection. 

The shocks located at Xi{t) ^^'^ Xr{t) must either be increasing or decreasing 
shocks. We will examine each of the possibilities and show that each leads to 
it(a;, t) not being a reparameterization of the initial conditions. 

Case 1 Assume that the shock at xii^) is a decreasing shock and the shock 
at Xr(i) is a decreasing shock. Then a(x/(i)~, t) > b{xi{t)~^,t) and b{xr{t)^ ,t) > 
c{Xr{t)~^ ,t). If b{xi{t)'^ , t) > 6(xr (i) ^7 i), then by the transitive property a(x;(t)~, i) > 
c(xr(i)~'"j t) and this violates the Rankine-Hugoniot condition, as was mentioned 
above, and u{x,t) is not a weak solution. If b{xi{t)^ ,t) < b{xr{t)^ ,t), then 
a(x;(i)~,^) > ^(Xi(^)^j^) < i'{Xr(t)~ ,t), shows u(x,t) is not a reparameteriza- 
tion of initial conditions. See figure |5] 



Figure 5: With two decreasing shocks, u{x,t) cannot be a reparameterization. 
The circles represent the points chosen to prove that u{x,t) cannot be a repa- 
rameterization of initial conditions. 
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Case 2 Assume that the shock at xi (^) is a decreasing shock and the shock at 
Xr{t) is an increasing shock. Then a{xi{t)^,t) > b{xi{t)~^,t) and b{xr{t)~,t) < 
c{xrit) + ,t). Let 62 = min {b{xi{t) + ,t),b{xAt)-,t)) , then a(x/(i)",<) > b2 < 
c{xr(t)~^ shows u(x,t) is not a reparameterization of initial conditions. See 
figure |6] 



Figure 6: With a decreasing and an increasing shock, u{x,t) cannot be a repa- 
rameterization. The circles represent the points chosen to prove that t) 
cannot be a reparameterization of initial conditions. 



Case 3 Assume that the shock at xi{t) is an increasing shock and that the 
shock at Xr{t) is an increasing shock. Since a{xi{ti)~ ,ti) — c{xr{ti)~ ,ti), for 
at least a short period of time after ti, the left value of b{x,t) will be greater 
than a{xi{t),t) and c{xi{t),t) and the right value of b{x,t) will be lower. By 
choosing the points b{xi{t)'^,t) > b{xr{t)~,t) < c{xr{t)^ ,t), u{x,t) is not a 
reparameterization of the initial conditions. See figure [7j 

Case 4 Assume that the shock at Xii^) is an increasing shock and that the 
shock at Xr {t) is a decreasing shock. This case will be divided into two subcases. 



U 




X 
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Figure 7: With an increasing shock and an increasing shock, u(x, t) cannot be a 
reparameterization. The circles represent the points chosen to prove that u{x, t) 
cannot be a reparameterization of initial conditions. 



The first is that for all x € {xi{t) , Xr{t)) that b{x,t) > a{xi{t)~,t) and b{x,t) > 
c{xr{t)~^ ,t). If this is the case, then u{x,t) is proven to not be a weak solution 
by Lemma [42j 

The second case is that there exists an xi G (Xi(^)j Xr{t)), such that b{xi, t) < 
a{xiit)~ jt) orb{xi,t) < c{xrit)~ ,t). Since x;(i) is an increasing shock and Xr(t) 
is a decreasing shock, and a(x;(ti)~, ii) = c{xr{ti)~ for at least a short 
period of time after ti, the left and right value of b{x,t) will be greater than 
aixi{t),t) andc(x/(i),t). Thus if < a{xi{t)- ,t) or6(a:i,t) < c{xrit)-,t), 

the points b{xi{t)^;t) > < b{xr{t)^ ,t) show that u{x,t) is not a repa- 

rameterization of initial conditions. See figure [Sj 

Thus if u(a;, t) is a weak solution of the inviscid Burgers equation and a 
reparameterization of initial conditions, it cannot engage in spontaneous shock 
formation. 
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Figure 8: With a decreasing and an increasing shock, u{x,t) cannot be a repa- 
rameterization. The circles represent the points chosen to prove that u(x, t) 
cannot be a reparameterization of initial conditions. 



4.3.2 Shock Splitting 

Assume that u{x,t) is the entropy solution to the inviscid Burgers equation up 
to time ti where an existing decreasing shock splits into two or more at point 
X*. Thus we say after time ti, 

{a{x,t) iix<xi{t), 
b{x,t) if xi{t) < X < Xr{t) (56) 
C{x,t) if Xr{t) < X, 

where a{x, t), b{x, t), and c{x, t) are weak solutions to the inviscid Burgers equa- 
tion and x/(i) and Xrit) give the position of the leftmost and rightmost disconti- 
nuities formed during the shock splitting. As there is assumed to be an already 
existing decreasing shock at time ti, a{xiiti)^ ,ti) > c(xr (<i)^, ^i). There may 



be more than two shocks formed as seen in Equation (4.1 ), but we just need to 
examine the leftmost and right most. 
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For u{x, t) to be a weak solution we note that several things must be true. 
The speed of xi{t) ^^'^ Xr(t) are dictated by the Rankine-Hugoniot jump con- 
ditions to be 

For there to be shock splitting, there must be some interval {ti,t2), where J^Xr > 
^Xi- Thus for some interval (ti,t2), o,{Xi) > c(Xr) ^■'^d thus b{xi{t)^,t) < 
b{xr{t)~ ,t). Assume that t e (ti,t2) for the remainder of the subsection. 

The shocks located at xi{t) and Xr{t) must either be increasing or decreasing 
shocks. We will examine each of the possibilities and show that each leads to 
t) not being a reparameterization of the initial conditions. 



Case 1 Assume that the shock at xi {t) is a decreasing shock and the shock at 
Xr{t) is a decreasing shock. Then a{xi{t)^,t) > 6(x;(t)+,t) and b{xr{t)^,t) > 
c{xr{ty,t). We know that b{xi{t) + ,t) < b{xr{t)',t), and thus a{xi{ty,t) > 
6(x;(t)^,i) < b{xr(t)^ ,t), shows u{x,t) is not a reparameterization of initial 
conditions. See figure [9j 



Case 2 Assume that the shock at xi (t) is a decreasing shock and the shock at 
Xr{t) is an increasing shock. Then a{xi{t)~,t) > b{xi{t)~^,t) and c{xr{t)'^,t) > 
b{xr{t)-,t) > b{xi{t) + ,t). Thus the points a(x;W",^) > h{xi{t) + ,t) < c(x.(t)+,t), 
shows u{x, t) is not a reparameterization of initial conditions. See figure 10 



Case 3 Assume that the shock at Xi(^) is an increasing shock and that the 
shock at Xr{t) is an increasing shock. Then a{xi{t)^,t) < b{xi{t)^,t) and 
c{xr(t)~^ ,t) > b(xrit)~ ,t). From the Rankine-Hugoniot jump conditions, for 
the interval (^1,^2), o-ixT) > ^(x^) and thus b{xi{t)'^,t) < b{xr{t)~ ,t). This is 
a contradiction, so u{x, t) is not a weak solution. 



Case 4 Assume that the shock at Xii^) is an increasing shock and that the 
shock at Xr{t) is a decreasing shock. This case will be divided into two subcases. 
The first is that for all x £ {xi{t)i Xr{t)) that b{x, t) > a{xi{t)^ ,t). If this is the 



case, then u{x,t) is proven to not be a weak solution by Lemma 4.2 



The second case is that there exists an Xi S (Xi(^): Xr(^)), such that b{xi,t) < 
a{xi{t)^ ,t). Since Xii^) is an increasing shock b{xi{t)'^,t) > a{xi{t),t) and thus 
6(xi(t)+,i) > b{xi,t). Since b{xi{t) + ,t) < b{xr{t)-,t), the points b{xi{t) + ,t) > 
< b{xr{t)~ ,t) show that u{x,t) is not a reparameterization of initial 
conditions. See figure [TT| 

Thus if t) is a weak solution of the inviscid Burgers equation and a 
reparameterization of initial conditions, it cannot engage in spontaneous shock 
formation or shock splitting. 
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Figure 9: With two decreasing shocks, u{x,t) cannot be a reparameterization. 
The circles represent the points chosen to prove that u(x, t) cannot be a repa- 
rameterization of initial conditions. 



4.3.3 Entropy violating solutions are not reparameterizations of ini- 
tial conditions 

From the previous sections the following lemma can be established. 

Lemma 4.3 Let u{x, t) be a weak solution of the inviscid Burgers equation 
where the initial conditions satisfy condition B. If u{x, t) is reparameterization 
of initial conditions then it is the entropy solution. 



Proof Clearly this follows from the results in Sections 4.2| 4.3.1 and 4.3.2 



4.4 Convergence to the Entropy Solution 



Based on Lemma 4.3 the following theorem regarding the CFB equations con- 
verging to the entropy solution can be established. 
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Figure 10: With a decreasing and an increasing shock, u{x,t) cannot be a 
reparameterization. The circles represent the points chosen to prove that u(x, t) 
cannot be a reparameterization of initial conditions. 



Theorem 4.4 The solutions u" of the CFB equations converge to the entropy 
solution of the inviscid Burgers equation for initial conditions satisfying condi- 
tion B. 

Proof It has already been established that converges to a weak solution of 
the inviscid Burgers equation, u, in L\^^. In the existence uniqueness proof, it 
was established that is a reparameterization of initial conditions. Clearly if 
every is a reparameterization, then its limit u will also be a reparameteri- 
zation of initial conditions. Since u is a reparameterization of initial conditions 



and a weak solution to the inviscid Burgers equation, by Lemma 4.3 u must be 
the entropy solution. 

Now that we have established that the solutions of the CFB equations con- 
verge to the entropy solution of the inviscid Burgers equation for initial condi- 
tions satisfying condition B. The following section deals with how to regain the 
entropy solution for discontinuous initial conditions and why we believe that 
this result hold true for more general cases. 
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Figure 11: With a decreasing and an increasing shock, u(a;, cannot be a repa- 
rameterization. The circles represent the points chosen to prove that u{x, t) 
cannot be a reparameterization of initial conditions. 



5 Extension into discontinuous initial conditions 

Section |4] proves that the CFB equations will converge to the entropy solution 
for a specific set of initial conditions. This section explains the intuitive reason- 
ing on why it is suspected that the CFB equations will converge to the entropy 
solution for any continuous initial conditions and why it will not for discontinu- 
ous initial conditions. It then shows how the equations can be changed slightly 
to incorporate discontinuous initial conditions. Begin by examining a commonly 
examined problem for the inviscid Burgers equation. 
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5.1 Example of entropic and non-entropic behavior for the 
inviscid Burgers equation 

Consider the initial conditions 



if X < 

1 if x > 0. 



(58) 



The method of characteristics does not provide the value of u in the wedge 
< X < t as seen in figure 12 The entropy solution fills this wedge with the 



function u{x,t) — | with characteristics fanning out from the original discon- 
tinuity as seen in figure |13^ . This creates a rarefaction wave and eliminates 
the discontinuity after time t = 0. A non-entropic solution will allow the dis- 
continuity to continue to exist. It will fill the wedge with new characteristics 
which continuously originate from the discontinuity as time progresses as seen 
in figure [T3| 3. Thus the non-entropic solution creates new 'information' as time 
progresses. 

This problem embodies the essential behavior of entropic and non-entropic 
solutions and provides the basis for our reasoning in the following subsections. 




Figure 12: The area < a; < t is not filled by characteristics 



5.2 Convergence to entropy solution for all continuous ini- 
tial conditions 

In section [4] it was proven that for initial conditions satisfying condition B, that 
the solutions to the CFB equations converge to the entropy solution. It is the 
conjecture of this paper that the solutions to the CFB equations converge to 
the entropy solution for all continuous initial conditions. As mentioned above 
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X X 

(a) (b) 

Figure 13: This is how entropic and non-entropic solutions fill the empty wedges, 
a) The entropic solution fills the wedge with a rarefaction wave, b) The non- 
entropic solution maintains the discontinuity by creating new characteristics. 



a non-entropic solution will create new characteristics, or new 'information' as 
time progresses. The solutions to the CFB equation do not. The existence 
and uniqueness theorem proven in our previous paper, :12J and restated here as 
Theorem |2.2[ established that the solution takes the form u{x,t) — UQ{(f>{x,t)), 
where (f){x, t) is an increasing function of x for any time, and 0) — x. This 
shows that no new information is being created in the CFB equations. Since the 
solutions to the CFB equations are converging to a weak solution to the inviscid 
Burgers equation and no new information is being created, it is reasonable to 
expect the solutions to converge to the entropy solution. 

5.3 Nonconvergence for discontinuous initial conditions 

Consider initial conditions that have an increasing discontinuity in them. The 
entropy solution to the inviscid Burgers equation creates a rarefaction wave from 
the discontinuity which takes on all the values spanned by the discontinuity. No 
new characteristics are formed, as all originate from the discontinuity at time 
t = Q, but u{x^t) now has values that did not originally exists in uq. As shown 
in the existence and uniqueness theorem for the CFB equations, the solutions 
to the CFB equations must have only the values found in wq. Thus for initial 
conditions containing an increasing discontinuity, the CFB equations will not 
converge to the entropy solution. An example of this can be found in our 
previous paper [12], in section 6, where a traveling wave solution to the CFB 
equations can be seen to converge to a non-entropic weak solution. For this 
reason, we eliminate discontinuous initial conditions from the admissible class 
of initial conditions. 

5.4 Conjecture 

Based on the reasoning in the previous two subsections, we present the following 
conjecture. 
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Conjecture 5.1 The solutions u" of the CFB equations converge to the entropy 
solution of the inviscid Burgers equation for continuous initial conditions as 

Assuming this conjecture is true, there is still the matter of discontinuous initial 
conditions. The following subsection creates a new system that if the conjecture 
is true will converge to the entropy solution for all bounded initial conditions. 

5.5 Regaining discontinuous initial conditions 

In regarding discontinuous initial conditions, begin by assuming that for all 
initial conditions the solutions to the CFB equations converge to the entropy 
solution. Then if the initial conditions limit to the discontinuous initial 
conditions in L^^^, at the same time as a ^ then the solutions will converge 
to the entropy solution for the discontinuous initial conditions. To prove this 
we use a theorem proven by Oleinik 15j presented here. 

Theorem 5.2 Let t) be the entropy solution for the inviscid Burgers equa- 
tion with initial conditions u"{x, 0) = Uq{x) and Uq{x) < m for all n. Let 



for n — > oo for any compactly supported continuous function f{x). Then the 
sequence u"'{x,t) converges for n ^ oo to the entropy solution u{x,t) in L\^^ 
with initial conditions w(x,0) = uq{x). 

This theorem is employed in proving the following theorem. 

Theorem 5.3 Let u"'" be solutions to the CFB equations with initial conditions 
■u"'"(a;,0) = Uo{x). Let Ug(x) converge to Uq{x) in as n ^ oo. Let 
be the entropy solution to the inviscid Burgers equation with initial conditions 
u'^{x^O) = uo{x). If for initial conditions uo{x) € C'^, the solutions to the 
CFB equations converge to the entropy solution in L\^^ , then u"'" converges to 
u'^{x, t) in L\^^ as n ^ oo and a — > 0. 

Proof Let be a compact subset of M x [0, T]. For u"'" to converge to u'^{x, t) 



Let u"'{x, t) be the entropy solution to the inviscid Burgers equation with initial 
conditions w"(a:,0) — Uq{x). We have assumed that 




in Ll 




(59) 




(60) 



From Theorem 15 . 21 we know that 




(61) 
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Thus employing the triangle inequality we find 

hm / / \u^'^ -u-\ < lim„^oo a^o I In I""'" - ""I 

+ lim„^oo a^o / In - "1 = (62) 

Using Theorem |5.3| it is easy to see that the solutions to the initial value 
problem 

ut + iu* g°')u^ = (63) 
u(x,0) = uo (64) 

will converge to the entropy solution of the inviscid Burgers equation with any 
initial condition uo{x) as a — > 0. This scheme can handle discontinuous initial 
conditions, providing a greater usefulness. 

6 Numerics 

Section |5] proposes that Equations ( 63 1 and ( 64 ) are a new system for the con- 
vectively filtered Burgers equation that is expected to converge to the entropy 
solution of the inviscid Burgers equation as a ^ for all bounded initial condi- 
tions. This section runs some numerical simulation of the proposed system and 
shows evidence of convergence to the entropy solution. 

6.1 The entropy solution 

The specific initial condition being examined is the indicator function for the 
interval (1,2) or 

Ji ifxe (1,2) 
I otherwise. 

For the entropy solution to the inviscid Burgers equation, the right side of the 
initial pulse will form the standard right traveling shock and the left side will 
form a rarefaction wave. At time t = 2, the rarefaction wave meets with the 
shock front the shock front begins to decrease in amplitude and speed. For time 
t < 2 the entropy solution for the given initial conditions is 

'0 ifa;<l 

K-,t) = { ^ if.-e(U + t) 

^ ' ^ 1 if X e {l + t,2 + M) ^ ' 

S) ifx>2+.5i. 

For time t > 2 the entropy solution is 

{0 if a; < 1 
^ if (l,(2t)5 + 1) (67) 
if a; > (2t)3 + 1. 

It is to this solution that the CFB equations' solutions are compared. 



30 



6.2 Description of Numerical Methods 

Holm and Staley performed successful simulations of the CFB equations with 
the Helmholtz filter, using a pseudospectral method [ST]. For this paper a very 



similar method is used. With the Helmholtz filter, Equations ( 63 1 and ( 64 ) can 
be written as 

d_ d 3 2/ 2d'V' d ^_ ,2 



-u + — — — = —-a 



?2(x,0) = (uo*5")*ff". (69) 
It is these equations that are numerically simulated. 



Equation (68 1 is advanced through time with an explicit, Runge-Kutta- 
Fehlberg predictor/corrector (RK45). The initial timestep is chosen low enough 
to achieve stability, and is then varied by the code using the formula 

h.+i=ihj '\ y. (70) 

\\\Ui-Ui\\2j 

Thus the new time step is chosen from the previous time step and the amount 
of error between the predicted velocity, u and the corrected velocity -u. The 
relative error tolerance was chosen at £ = 10^^ and the safety factor 7 = 0.9. 

Spatial derivatives and the inversion of the Helmholtz operator were com- 
puted in the Fourier domain. The velocity was converted into the Fourier do- 
main using a Fast Fourier Transform, multiplied by the appropriate term and 
then converted back into the physical domain. This pseudospectral method of 
calculating the derivative was chosen to reduce artificial viscosity. 

In Holm and Staley's method spatial derivatives were conducted using a 
fourth-order finite difference and an artificial viscosity was applied to the high 
wave modes to prevent aliasing errors [31j . Because the simulations are ad- 
dressing convergence to the entropy solution, as little artificial and numerical 
viscosity as possible is desired. For this reason derivatives were done in the 
Fourier domain and no artificial viscosity was introduced. 

The simulations were done at the resolution of 2^^ = 65536 grid points. 
Aliasing errors occurred, but did not introduce significant amounts of error in the 



short time the simulations were run. Figure 14 shows the spectral energy of the 
simulation for a = 0.02 at various times. This was the worst case of aliasing error 
and it can be seen that the error does not reach more than approximately 10"^^, 
at t = 3 which is approximately 100,000 timesteps. It should be noted that 
simulations using Holm and Staley's artificial viscosity prevented this aliasing 
error, with little noticeable effect on the solution. Simulations with the artificial 
viscosity have been conducted and produce the same general results presented 
in the following sections. 

6.3 Results 

Nine different simulations were conducted with a = 0.02, 0.03, 0.10. The CFB 
equations showed behavior mirroring that of the entropy solution. A traveling 
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,qZ5| 1 , 1 1 1 ,g2s| 1 , 1 1 1 

10° 10' 1 10' 10' 10' 10° 10' 1 10' 10' 10' 

a it a it 

(c) (d) 

Figure 14: The error caused by aliasing is visible in the lower part of the energy 
spectrum. These are snapshots of the spectral energy at times t = 0.5,1,2,3. 
The aliasing error is seen to be propagating up into lower wave modes, but at 
time t = 3, the error is roughly capped off by 10~^^. The dashed lines give a 
reference -2 and -6 slope. The spectral energy slope changes at approximately 
^ as is expected. 



shock front and a rarefaction wave was seen. Figure [15] compares the CFB 
simulations for a — 0.02 to the entropy solution at times t = 0,1, 2, 3. In figure 
|15^ the difference in initial conditions can be seen with the entropy solution 
beginning with discontinuities and the CFB simulation having smoothed initial 
conditions. 

To evaluate the convergence of the CFB equations' solutions to the entropy 
solution the Li norm of the error between the CFB equations' solution and 
the entropy solution was taken. Figure [16] plots a versus the error at times 
t = 0,1,2,3. At each time the error appears to be approaching zero linearly. 
Thus numerical evidence suggests that the Equations ( 63 1 and ( 64 1 will converge 
to the entropy solution of the inviscid Burgers equation for initial conditions with 
discontinuities. 
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(c) (d) 

Figure 15: This figure compares the entropy solution with the solution to the 
CFB equations for a — 0.02. It is easy to see that the CFB equations' solution 
is capturing both the rarefaction wave and the shock front behavior. 



7 Conclusion 

Conservation laws can often have multiple weak solutions of which there is one 
physically relevant solution, known as the entropy solution. It is important that 
any regularization of these conservation laws reflect the physical phenomenon 
they are meant to address. Thus it is important that the solutions to such reg- 
ularizations converge to the entropy solution. The convectively filtered Burgers 
equation has been shown to regularize the inviscid Burgers equation. This paper 
now shows that for a certain class of initial conditions this regularization will 
converge to the entropy solution. It has also provided a method for extending 
this convergence to a large class of initial conditions including discontinuities. 
These results are a crucial step in extending the use of the convectively filtered 
method into popular use and perhaps an extension into the Euler equations. 
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Figure 16: The Li norm of the error between the CFB equations' solution and 
the entropy solution. The error is displayed for four different values of t. t = 

, t = 1 , t = 2 , and t — 3 . The error approaches zero 

roughly linearly as a 0. 
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